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1 Introduction 



In the Standard Model (SM), the parameter 



P 



w 



\ cos^ Qw 



m 



where m^y and mz are the masses of the and gauge bosons, respectively, and 9w 
is the weak mixing angle, gives the relative strength of the neutral-current and charged- 
current interactions in four-fermion processes at zero momentum transfer p]. At tree 
level p is equal to one, and it remains one even if additional scalar SU{2) doublets, with 
hypercharge ±1/2, are added to the Sm]^ At one-loop level, the vacuum-polarization 
effects, which are sensitive to any field that couples either to the or to the 
produce the vacuum-polarization tensors {V = W, Z) 



(2) 



where is the four-momentum of the gauge boson. Then, deviations of p from unity 
arise, which are determined by the self-energy difference [H [2] 
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The precise measurement [3], at LEP, of the and Z° self-energies is in striking agree- 
ment with the SM predictions [1] and provides a strong constraint on extended electroweak 
models. For instance, one can constrain the two-Higgs-doublet model (2HDM) in this 
way [HIE]. 

In this paper we are interested in the contributions to the p parameter generated by 
an extension of the SM. Therefore, we define a Ap which refers to the non-SM part of the 
quantity ([3]): 
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SM extension 
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The SM contributions to the quantity ([3]) are known up to the leading terms at three-loop 
level [7]. However, the consistent SM subtraction in equation (jlj) only requires the one- 



loop SM result. In the same vein, we are allowed to make the replacement 
in equation (jlj), writing instead 
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SM 



Here and in the following, we use the abbreviations cw = cos Ow^ sw = sin 9w- 

At one loop, the contributions of new physics to the self-energies constitute intrinsically 
divergent Feynman diagrams, but the divergent parts cancel out among different diagrams. 



^Other scalar SU{2) x U{1) representations are also allowed, as long as they have vanishing vacuum 
expectation values. 
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between A^w (0) and c^Azz (0), and also, eventually, through the subtraction of the SM 
contributions laid out in equation ([5]). If the new-physics model is renormalizable, then 
Ap is finite. The cancellations finally leave either a quadratic or a logarithmic dependence 
of Ap on the masses of the new-physics particles. The pronounced effects of large masses 
is what renders the parameter Ap so interesting for probing physics beyond the Standard 
Model. 

The functions Ayv {(f) contain more information about new physics than the one just 
provided by Ap. In fact, for new physics much above the electroweak scale, a detailed 
analysis of the so-called "oblique corrections" lead to the identification of three relevant 
observables, which were called 5, T and f/ in [8] and ei, €2 and £3 in [9J@ While these two 
sets of observables differ in their precise definitions, the quantity of interest in this paper 
is simply 

Ap = aT = ei, (6) 

where a = e^/ (47r) = g^s"^/ (in) is the fine-structure constant. 

It is not straightforward to obtain a bound on Ap from electroweak precision data. 
One possibility is to add the oblique parameters to the SM parameter set and perform 
fits to the data. However, since the SM Higgs-boson loops themselves resemble oblique 
effects, one cannot determine the SM Higgs-boson mass rrih simultaneously with S and 
T To get a feeling for the order of magnitude allowed for Ap, we quote the number 

T = -0.03 ±0.09 (+0.09), (7) 

which was obtained in [4J by fixing U = 0. For the mean value of T, the Higgs-boson 
mass rrih = 117 GeV was assumed; the mean value in parentheses is for rrih = 300 GeV. 
Equation (CD translates into Ap = -0.0002 ± 0.0007 (+0.0007). 

There is a vast literature on the 2HDM — see [llj for a review, [I2] for the renormal- 
ization of the model, [TSj [H] for the possibility of having a light pseudoscalar compatible 
with all experimental constraints, and |15i, 116j , and the references therein, for other var- 
ious recent works. However, just as the 2HDM may differ significantly from the SM, a 
general multi-Higgs-doublet model may be quite different from its minimal version with 
only two Higgs doublets [H]. Three or more Higgs doublets frequently appear in models 
with family symmetries through which one wants to explain various features of the fermion 
masses and mixings; for some examples in the lepton sector see the reviews in |18j . 

In this paper we present a calculation of Ap in an extension of the SM with an arbitrary 
number of Higgs doublets and also, in addition, arbitrary numbers of neutral and charged 
scalar SU{2) singlets. Our results can be used to check the compatibility of the scalar 
sector of multi-Higgs models with the constraints resulting from the electroweak precision 
experiments. 

Recently, there has been some interest in "dark" scalars [191 [20]. These are scalars that 
have no Yukawa couplings, and are thus decoupled from ordinary matter. Furthermore, 
they have no vacuum expectation values (VEVs) and therefore display truncated couplings 
to the gauge bosons. However, they would have quadrilinear vector-vector-scalar-scalar 

^For new physics at a mass scale comparable to the electroweak scale three more such "oblique pa- 
rameters" have been identified in [lOj . 
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and trilinear vector-scalar-scalar (but no vector-vector-scalar) couplings, and would thus 
also contribute to, and be constrained by, Ap. 

The plan of the paper is as follows. In Section 2 we present a description of our 
extension of the SM and the final result of the calculation of Ap; this section is self- 
consistent and the result can be used without need to consult the rest of the paper. The 
details of the calculation are laid out in Section 3. The application of our Ap formula to 
the general 2HDM is given in Section 4. The summary of our study is found in Section 5. 



2 The model and the result for Ap 
2.1 The model 

We consider an SU{2) x U{1) electroweak model in which the scalar sector includes rid 
SU{2) doublets with hypercharge 1/2^ 



k = 1,2,..., rid. 



Moreover, we allow the model to include an arbitrary number and variety of SU (2)-singlet 
scalars; in particular, complex SU{2) singlets with hypercharge 1, 

X/, J = l,2,...,n, (9) 

and Un real SU{2) singlets with hypercharge 0, 

Xl /= l,2,...,n„. (10) 

In general, our model may include other scalar fields, singlet under the gauge SU{2), with 
different electric charges. 

The neutral fields are allowed to have vacuum expectation values (VEVs). Thus, 

<0K|0> = ^, (11) 
(0|x?|0> = ui, (12) 
the Vk being in general complex. (The ui are real since the are real fields.) We define as 

1 /2 

usual V = (^^ii \vkf) ' ~ 246 GeV. Then, the masses of the and gauge bosons 
are, at tree level, mw = gv/2 and niz = mw/cw, respectively|3 We expand the neutral 
fields around their VEVs, 

= -^{vk + d% (13) 



■^Equivalently, we may consider the model to contain SU(2) doublets with hypercharge —1/2, since 

(f>k = iT2(j)l = 



is also a doublet of SU{2). 

'''Since the neutral singlet fields carry no hypercharge, their VEVs ui do not contribute to the masses 
of the gauge bosons. 
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X'l = u, + xr- (14) 

Altogether, there are n = + n^. complex scalar fields with electric charge 1 and 
m = 2nrf + Tin real scalar fields with electric charge 0. The mass matrices of all these 
scalar fields will in general lead to their mixing. The physical (mass-eigenstate) charged 
and neutral scalar fields will be called S'^ (a = 1, 2, . . . , tt.) and {h = 1,2, ... ,m), 
respectively. Note that the fields are real. We use rria to denote the mass of and 
fj,b to denote the mass of S^. We have 

n 

= J2^kaS^, (15) 
a=l 
n 

Xt = Y^^^aSt, (16) 
a=l 
m 

Y.'^kA (17) 



b=l 



X?' = (18) 
b=i 

the matrices U , T , V and R having dimensions x n, ric x n, rid x m and n„ x m, 
respectively. The matrix R is real, the other three are complex. The matrix 

is 77. X n unitary; it is the matrix which diagonalizes the (Hermitian) mass matrix of the 
charged scalars. The real matrix 



V= ImV (20) 




is m X m orthogonal; it diagonalizes the (symmetric) mass matrix of the real components 



of the neutral-scalar fieldsl 

There are in the spontaneously broken SU (2) x U (1) theory three unphysical Goldstone 
bosons, G"^ and G^. For definiteness we assign to them the indices a = 1 and 6=1, 
respectively: 

St = G^, (21) 
5° = G°. (22) 

Thus, only the with a > 2 are physical and, similarly, only the with b >2 correspond 
to true particles. In the general 't Hooft gauge that we shall use in our computation, the 
masses of G^ and G^ are arbitrary and unphysical, and they cannot appear in the final 
result for Ap. 



■^Our treatment of the mixing of scalars is inspired by 
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2.2 The result 



As we shall demonstrate in the next section, the value of Ap in the model outlined above 
is 

^ L a=2 b=2 

m— 1 m 

-EE [lm(V^V)jV(p^,p^,) (23b) 



b=2 b'=b+l 
n— 1 n 



2E E \iU'U)jF{mlml,) (23c) 



=a+l 



+3 5^ [Im [VW] ^ J ^ [F (m|, p^) - F (m^^, p^)] (23d) 

6=2 

-3 [F (m|, m^) - F (m^, ml)] } , (23e) 

where rria, ma' denote the masses of the charged scalars and p^, pf,/ denote the masses 
of the neutral scalars. The term (123bp contains a sum over all pairs of different physical 
neutral scalar particles and S^, ; similarly, the term fl23cp contains a sum over all pairs 
of different charged scalars, excluding the Goldstone bosons G^, i.e. 2 < a < a' <n. The 
term fl23el) consists of the subtraction, from the rest of Ap, of the SM result — mh is the 
mass of the sole SM physical neutral scalar, the so-called Higgs particle. 

In equation fl23l) . the function F of two non- negative arguments x and y is 

{x + y xy X 
m — <^= a; ^ y, 
2 x-y y (24) 
<^ X = y. 

This is a non-negative function, symmetrical under the interchange of its two arguments, 
and vanishing if and only if those two arguments are equal. This function has the impor- 
tant property that it grows linearly with ma.x{x,y), i.e. quadratically with the heaviest- 
scalar mass, when that mass becomes very large. Unless there are cancellations, this leads 
to a quadratic divergence of Ap for very heavy scalars (Higgs bosons). 

If there are in the model any S'f/(2)-singlet scalars with electric charge other than 
or ±1, then the existence of those scalars does not contribute to Ap, they do not modify 
equation (l23!l . at one- loop level, in any way. 

A simplification occurs when there are in the model no SU (2)-singlet charged scalars 
X^- In that case, there is no matrix T, hence the matrix U is unitary by itself, and the 
term fl23cl) vanishes. 

When there are in the model no St/ (2) -singlet neutral scalars there is no matrix 
R, hence Re (^^V^)^^, = {ReV^ ReV + ImV^ 1^^)^^^, = 6w Then, in the terms fl23bl) 

and fl23dp one may write | (^^^)^^/|^ instead of [im (^^^)^^/] 



Thus, in an ra^f-Higgs-doublet model without any scalar singlets, one has simply 

2 ( "-d 27id 



w 



, a=2 6=2 
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271^ — 1 2nd 

6=2 fe'=6+l 
2nd 

6=2 

-3[F {ml,ml) - F {m'^,ml)]} . (25) 

Our general results have been checked to be consistent with specific results for Ap in a 
few models. These include the results for both the CP conserving version [5l [13] and the 
CP non-conserving version [16] of the 2HDMj!l It has also been checked against a model 
containing one doublet and one scalar singlet |23j . 



3 Derivation of the result 

This section contains the derivation of equation (|23l) . It may be skipped by those who 
are not interested in the details of that derivation. 



3.1 The Lagrangian 

We use the conventions of [22j. The covariant derivative of the doublets is 



2c 



w 



(26) 



and the covariant derivative of the charged singlets is 

2 



(27) 



The covariant derivative of the neutral singlets is, of course, just identical with their 
ordinary derivative. We use the unitarity of U in equation (JT9l) . in particular 



{2i 



We also use the orthogonality of V in equation (I2U|) to arrive at the gauge-kinetic La- 
grangian 

nd nc ^ n„ 

{D^<P,y {DM + [D^x]) {D,xt) + 7^Y1 i^'^i) i^^^i) 



k=l 



1=1 



n ^ m 



(29a) 



a=l 



6=1 



^There is some discrepancy between our result and the one presented in Section 4 of [TT]. 
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z^z 

+ml,W^'-W+ + m\ (29b 

n 

+im^ J2 [W~ {u^^U)J^S^ - {U^u)J^S^] (29c 

a=l 
m 

+mzZ^ J2 Im {oo^V) ^ d^S^, (29d 

6=1 

n 

- {emwA^ + gslmzZ^) [{uJ^U)^W-S: + {U^u)^ W^^S'] (29e 

a=l 

n 

+zeA, J2 {S^d^S- - S-d^S:) (29f 



a=l 

n 



[24'5aa' - {U^U)J {S^d^S-, - S;,d^S^) (29g 

a,a'=l 
m— 1 m 

n m 
a=l b=l 

+ (V^U)^^ W- {S'.d^S^ - S^d^S',)] (29i 
+, L^W^W^- + ^ ^) ^ 5° Re (^t^)^ (29j 

/ 2 2 \ n m 

+ (^ + ^ ^) E (^-^f).. sjsr (291 

2 " 

Yl iU'U)^,^S;,S^ (29m 



a,a '= 1 



+2e^^j2SaS: (29n 



2 

a=l 

n 

eg 



a,a'=l 



+ ^ ^ E [44^..' + (1 - 4.^H^) (t/^f/),J S~,St. (29p 



a,a'=l 



In lines fl29cl) - fl29ep and ([29]]) we have used an n^- vector uj defined by uj^ = Vk/v. By 
identifying lines (I29cp and (I29dp with the usual terms [22] mixing the and gauge 
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bosons with the and Goldstone bosons, respectively, 

imw {W-d^'G+ - W^d^'G-) + mzZ^d>'G^, 
we conclude that the components of the Goldstone bosons are given by [21 



Uki = — , hence T,i = 0, (30) 

V 

Vki = i hence Rn = 0. (31) 

V 

Therefore, we may rewrite line f l29ep as 

- {emwA^ + gsl.mzZ'') {W-G+ + W^G') (32) 

and line p9jP as 



- , (m^W^W^- + ^ ^) 5° Im {V^V) . (33) 

V Cw ^ / 

The sum starts at 6 = 2 because Im (yl"!^)^,^ = 0. 

If there are in the theory any S'?7(2)-singlet scalars S^'^ with electric charge ±Q other 
than or ±1, then those scalars do not mix with components of the doublets. Their 
covariant derivative is 



= + t + ieQA^S+'^. (34) 



This yields, in particular, the following two interaction terms in the Lagrangian: 

£=...+, - (35a) 

Cw 
\ Cw J 

3.2 The Feynman diagrams 

In our model, in the computation of the vacuum polarizations of the gauge bosons 
and there are four types of Feynman diagrams involving scalar fields: 

Type (a) diagrams: A scalar branches off from the gauge-boson line and loops back to the 
same point in that gauge-boson line — see figure [T](a). When the scalar is neutral, 
the relevant interaction terms in the Lagrangian are the ones in line (I29ip . for b' = b; 
but then the contribution to Ap vanishes, since one obtains H^^^y = c^^U^^z- When 
the scalar is charged, the relevant terms in the Lagrangian are those in line fl29ml) 
for ^'^lY and line p9pP for II^'^, in both cases for a' = a. 

Type (b ) diagrams: The gauge-boson line splits into two scalar lines which later reunite to 
form a new gauge-boson line — see figure [U^b). The relevant terms in the Lagrangian 



are those in line (129ip for H'^yy, and those in lines p9gP and (]29hp for II^'^- 
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(a) (b) (c) 



Figure 1: Three types of Feynman diagrams occurring in the calculation of the vacuum 
polarizations. 



Figure 2: Tadpole diagrams which do not contribute to Ap. 

Type (c) diagrams: A neutral scalar branches off from the gauge-boson line and loops to 
a later point in that gauge-boson line — see figure [T](c). The interaction terms in the 
Lagrangian responsible for these Feynman diagrams are those in expression (!33l) . 

Type (d) diagrams: A neutral scalar branches off, with zero momentum, from the gauge- 
boson line, and produces a loop of some stuff — see figure El These "tadpole" Feyn- 
man diagrams originate from the interaction terms in expression fl33l) . They yield 
a vanishing contribution to Ap since one obtains ^'^y^ = c^Ii^zz- Hence we may 
omit the tadpole diagrams altogether. 

3.3 Computation of the loop diagrams 

We use dimensional regularization in the computation of the Feynman diagrams. The 
dimension of space-time is d. An unphysical mass p is used to keep the dimension of each 
integral unchanged when d varies. We define the divergent quantity 

2 

div = ^— ^ - 7 + 1 + In (47rp^) , 

where 7 is Euler's constant. In the computation of type (a) Feynman diagrams the 
relevant momentum integral is 

= ^ ^ (div - In A) , (36) 

{2T^f k^-A + ie IQtt^ ^ ^ ' 

where A is the mass squared of the scalar particle in the loop. In order to compute the 
type (b) and type (c) Feynman diagrams we need first to introduce a Feynman parameter 
x, which is later integrated over from x = to x = 1. For type (b) diagrams we have 

P / ;j / dx ^ = \A (div - In A) 

J {2Tif Jo [k^-Ax-B{l-x)+ief 167r2^ ^ ' 
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+B (div -\nB) + F {A, B)] , (37) 



where A and B are the masses squared of the scalars in the loop, and the four-momentum 
q'^ of the external gauge-boson line is taken to obey = 0. Notice the presence of terms 
of the form A (div — In A) in both diagrams of types (a) and (b); we shall soon see that 
those terms cancel out in the computation of Ap, leaving only the F functions from the 
type (b) diagrams. For type (c) diagrams the relevant integral is 



/i 



A-d 



dx ■ 



ig'''' 1 



[A;2 -Ax-B{l-x)+ ie] 



167r2 A 
+FiA,B)]. 



A (div - In A) 



A + B 



(St 



This integral is symmetric under the interchange of A and B; equation fl38l) presents a 
seemingly asymmetric form, but it is in fact symmetric. The reason for expressing the 
integral in this way is that, due to cancellations, only the terms F {A, B) survive in the 



end. 



3.4 The contributions to A/? from diagrams of types (a) and (b) 



Using fl29ml) and we see that the contribution to Ay/w {q^) of type (a) Feynman 
diagrams with charged scalars in the loop is 



^WW [1 ) - 
In the same way, using ( |29pP , 



327r2 



Y.iU'U)^ml (div-lnm^). 



(39) 



a=l 



4t (^^) 



327r2c2, 



[4/^ + (1 - Asl) (f/tf/) J ml (div - Inm^) . (40) 



a=l 



Proceeding to the type (b) Feynman diagrams, from fl29ip and (P7|) we find that 



9' 



647r2 



a=l 6=1 

+pI (div - In pI) + F [ml, pi)] 

'2Y,{U'U)^^ml{di^-\nmD 



a=l 



We have used 



647r2 

m 

6=1 

n m 

a=l 6=1 

n 

a=l 



bb ' 



(41a) 
(41b) 
(41c) 

(42) 
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which follows from the unitarity of U, i.e. from [21] 
We have also used 

m 
b=l 

which follows from the orthogonality of V, i.e. from |2T] 

ReVReV^ = ImV ImV^ = l„,x„„ 
ReVlmV^ = ImV ReV^ = On.xn,- 

Considering now the self-energy of the boson, we find 

-,2 



^ a,a'=l 



X [ml (div - Inm^) + ml (div - Inm^,) + F {ml, ml)] 

m— 1 m 
6=1 b'=b+l 

X [fJ'l (div - Inul) + nl (div - lufil) + F {filfil)] } 

,2 



n— 1 n 
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2E E m'^)J F«<<') 

a=l a'=a+l 

+2 5^ [4.^^ + (1 - Asl) (f/tf/) J (div - Inm^) 

a=l 
m— 1 m 

+ E E 



b=l 6'=6+l 
m 



+ E(^^^)../^ndiv-ln/.^) 



b=l 



We have used 



b'=l 



which follows from equations fH5l) . 

Putting everything together, we see that 

the A^^^y (g^) of equation fl39|) cancels out the line (14 lap of A^y^ (0); 
the A^z {(f) of equation (gO]) cancels out the line fl46bl) of (0); 



the line fl41bl) of ^4^^^^^ (0) cancels out the line fl46dl) of (0) the subtraction A 
CwAzz- 
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In this way we finally obtain 



<t^^ (0) - c^4T (0) = i^AJ:J:\iU'V)JF{ml^^l) (48a) 

L a=l 6=1 

n—1 n 

E \iU'U)jF{mlml,) (48b) 

a=l a'=a+l 
m—1 m 



6=1 b'=b+l 



The positive term fl48al) originates from A\^yy while the negative terms (]48b|) and fl48cl) 
come from A^^- 

If there are in the electroweak theory any scalar SU{2) singlets with electric charges 
other than or ±1, then the relevant terms in the Lagrangian are those in equation ( |35l) . 
The term fl35bl) generates a type (a) Feynman diagram which exactly cancels the type (b) 
Feynman diagram generated by the term fl35al) FI Thus, scalar SU (2) singlets with electric 
charge different from and ±1 do not affect Ap at all. 

The sums in equation ( l48l) include contributions from the Goldstone bosons = Sf 
and = 5*°. These Goldstone bosons have unphysical masses mi and pi, respectively, 
which are arbitrary in a 't Hooft gauge. The terms which depend on those masses are, 
explicitly, 

\{U^V)JF{ml,^,l) (49a) 

m 

+Ei(f^^^)i.r^«'"') (49b) 

6=2 

n 

+EI(f^'^)air^K/^D (49c) 

a=2 

n 

-2j2\{U'U)jF{mlml) (49d) 

a=2 
m 

-5^[lm(l^V)jV(;.?,;.^). (49e) 

6=2 

One may eliminate some of these terms by using equations (!30|) and (!3T|) . Indeed, 
(^WU)^^ = — (T^T)_^^ = when a 7^ 1, because Tji = for any j; also, (^WV)^^ = 
i (f/^f/)^_^ = for a 7^ 1. Therefore, the terms fl49cp and (]49dp vanish. In the term (149 al) . 
{WV)^^ = i. In the term fl49bD one may write 

(t/ty)^^ = z {V^V)^^ = -Im {V^V)^^ ^ M 1, (50) 

since Re {V^V)^^ = {ReV^ RelZ + ImV^^ ImV)^^ = - {R^R)^^ = 0. In this way, the 



^This cancellation is analogous to the one between equation pp]) and line (|46bp . 
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terms (H^ are reduced to 



F (ml^^l) + Y: [Im (yty) [F [ml ^^l) - F (/.?,/.^)] . (51) 

6=2 

The term F {ml, fif) is independent of the number of scalar doublets and singlets, hence 
it is eliminated when one subtracts the SM result from the Multi-Higgs-doublet-model 
one. The other terms in the expression ( ISTI) are cancelled out by the diagrams of type (c), 
as we shall see next. 



3.5 The contributions to Ap from diagrams of type (c) 

To compensate for the unphysical masses of the Goldstone bosons, the propagators of 
gauge bosons and with four-momentum k'^ are, in a 't Hooft gauge, 

^^^^ ' +(-9.^ + ^]t^^, (52) 



k"^ — mf \ J k"^ — ' 

k^ky i I k^kii^ i /"cqn 



m| k'^ — fil \ m| / k"^ — m 



respectively, i.e. they contain a piece with a pole on the unphysical masses squared 
and /if, respectively. 

Using these propagators to compute the type (c) Feynman diagrams, one obtains 

2 

^ww (0) = ^ E (^'^) 16] ' \r< (di^ - 1^ <) - ^^w (div - In m^) 

6=2 

+2 (m^, + ^^l)-F (m?, /i^) - 3F (m^^, , (54) 

2 ™ 

^^'^ = ^ufv^ E (^'^) 16] ' [-/^? (di^ - In /^D - 3m| (div - In m|) 



6=2 



+2 (m| + /^^) - F (/x^, _ 3^ (^1^ ^2^] _ ^55) 

The factors 3 originate in a partial cancellation between the contributions from the pieces 
—g^u and k^k^/my in the propagator of the gauge boson V, the former contribution being 
four times larger than, and with opposite sign relative to, the latter one, cf. equations (1571) 
and fl55]) . Performing the subtraction relevant for Ap, one obtains 

2 ™ 



A'^wi0)-cUzU0) = ^EP-(^'^)i6]- 



6=2 

X [—ml (div — + ^\ (div — In /if) (56a) 

— 3m^ (div — Inm^) + 3m| (div — lnm|) (56b) 

+2 {ml, - ml) (56c) 

-F(mf,/i^)+F(/if,/ia (56d) 

-3F (m^, /i^) + 3F (m|, /i^)] . (56e) 
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The terms fl56ap - fl56cp are independent of the number of scalar doublets. They disappear 
when one subtracts the Standard-Model result from the multi-Higgs-doublet-model one, 
since 

m 
6=2 

The terms ( ]56dl) . which involve the masses of the Goldstone bosons, cancel out the terms 
in ( ISTl) . except the first one, which is cancelled by the subtraction of the SM result. 
We have thus finished the derivation of equation fl23|) for Ap. 



4 The 2HDM and the Zee model 



In this section we give, as examples of the application of our general formulae, the expres- 
sions for Ap in the 2HDM and also in the model of Zee [23] for the radiative generation 
of neutrino masses, which has one singly charged SU{2) singlet together with the two 
doublets. 

In the study of the 2HDM it is convenient to use the so-called "Higgs basis," in which 
only the first Higgs doublet has a vacuum expectation value. In that basis. 



H + tG^) 



(58) 



Here, = and G^ = are the Goldstone bosons, while 5*^ is the physical charged 
scalar, which has mass 1712. Thus, the matrix U, which connects the charged components 
of 01 and 02 to the eigenstates of mass, is in the Higgs basis of the 2HDM equal to the unit 
matrix. On the other hand, H, R and J, which are real fields, must be rotated through 
a 3 X 3 orthogonal matrix O to obtain the three physical neutral fields 5*2 3 



i,3,4- 







(SI 


R 




SI 


• ) 




\ SI 



(59) 



Without lack of generality we choose det O = +1. Thus, the 2x4 matrix V, defined 
through 

/ GO \ 

H + iG° 
R + il 



V 



SI 



SI 
\SIJ 



(60) 



IS 



V 



o 



11 



o 



12 



o 



13 



021+i031 022 +iOz2 023 + zO: 



^33 



Therefore, 



/ 


1 


-^0n 


-iOi2 




\ 






1 




-iOi2 






iOi2 




1 






V 


^Oi3 


lOl2 




1 


/ 



(61) 



(62) 
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The value of Ap in the 2HDM is therefore, using our formula in equation (1251) . 



9^ 



647r2 



w 



J2 - Ol_^) F [ml fil) 



6=2 



+3 Ol,-i [F (m|, _ F {ml., ^l) - F (m|, ml) + F [m^, ml)] I , (63) 

b=2 ) 

where /i2,3,4 denote the the masses of 5*2 3 4, respectively, while m^ is the mass of the Higgs 
boson of the SM. Equation (IMj) reproduces, in a somewhat simplified form, the result for 
Ap in the 2HDM previously given in |16j . 

A special case of the 2HDM is the model with one "dark" scalar doublet. This means 
that a second doublet is added to the SM, but that doublet has no VEV and it does not 
mix with the standard Higgs doublet pi9]. We should then identify H with the usual 
Higgs particle. Thus, On = 1 and p2 = ^h- Equation fp5]) then simplifies to 



Ap 



J2F{ml,fil)-F{fil,fil) 



b=3 



(64) 



This quantity is small if the three masses m2, ps and p4 are close together. Notice that in 
this case of a "dark" scalar doublet there are no vector-vector-scalar couplings involving 
the additional doublet, hence Ap stems exclusively from type (a) and type (b) Feynman 
diagrams. 

In the model of Zee there is, besides the two scalar SU (2) doublets 



{v + H + zG") /V2 



H+ 
(R + iI) /V2 



(65) 



also one scalar SU(2) singlet x"*" with unit electric charge. Therefore there is a 2 x 2 
unitary matrix K such that 



K 



St 



(66) 



where S2 and are the physical charged scalars, which have masses m2 and ms, respec- 
tively. So, now the matrix U of equation (IT^ is 



U 



1 

i^n ^12 



(67) 



so that 



\Kii\ 
KuK 



ll-f^l2 







(6^ 
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Equations (15T|) and fl^ retain their validity, and 

/ t On Ou Oi3 \ 

U^V=\ Kl,i02l+^0sl) Kl,{022 + ^0^2) K*,, {O23 + lO^s) . (69) 
\0 Kl^{02l+^0sl) Kl^{022 + ^0^2) K*,^ {O23 + lO^s) J 

Therefore, using our general formula (|23|1 for Ap, we see that, in the model of Zee, 

W I b=2 a=2 

-2\KuK,2fF{mlml) 

+3 O'b-i [F K, fil) - F {m%,, fil) - F (m|, m^) + F (m^^, m^)] i . (70) 

6=2 J 

5 Summary 

In this paper we have derived the formula for the parameter Ap, as defined in equation (jlj), 
in an extension of the Standard Model characterized by an arbitrary number of scalar 
SU{2) doublets (with hypercharge ±1/2) and singlets (with arbitrary hypercharges) . Our 
formalism is completely general, using only the masses of the scalars and their mixing 
matrices, which ensures that our formulae are always applicable. The computation has 
been carried out in a general gauge, thereby demonstrating that the final result is 
independent of the masses of the unphysical scalars. We have also explicitly demonstrated 
that all infinities cancel out in the final result for Ap. In order to ease the consultation of 
this paper, the formulae for Ap given in Section 2 have been completely separated from 
their derivation presented in Section 3. Our results can be applied either to check the 
viability of a model or to constrain its parameter space, by comparing the Ap, calculated 
in that model, with numerical bounds on Ap obtained from a fit to precision data — 
for instance, the bound ([7]) found in [1]. As an illustration of our general formulae, in 
Section 4 we have worked out the specific cases of the two-Higgs-doublet model, with and 
without one extra charged scalar singlet. 
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